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Abstract. We present results for the wave functions and the screening mass spectrum for quantum numbers 
++ , 1 and 2 ++ in the three-dimensional 5 , (7(2)-Higgs model near to the phase transition line below 
the endpoint and in the crossover region. Varying the 3-D gauge couplings we study the behaviour along a 
line of constant physics towards the continuum limit in both phases. In the crossover region the changing 
spectrum of screening states versus temperature is examined showing the aftermath of the phase transition 
at lower Higgs mass. Different to smearing concepts we used large sets of operators with various extensions 
allowing to identify wave functions in position space. 
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1 Introduction 

During the last years, due to efforts of three groups Jl|-§ 
using the 31? approach, various aspects of the high temper- 
ature electroweak phase transition have been explored in 
the 5'C/(2)-Higgs model, with increasing Higgs mass Mh- 
The interface tension has been seen falling over three or- 
ders of magnitude (Htl expressed as a/Tj?), whereas 
the latent heat goes to zero roughly linearly with Mh- 
The continuum characteristics of the transition is difficult 
to obtain with the same precision by 4-D Monte Carlo sim- 
ulations 0, although the results are consistent with each 
other where they can be compared ^,||. Approaching the 
W mass Myy, the character of the transition is now known 
to change drastically. It is the spectrum of physical excita- 
tions in this parameter region that we are going to study 
here. 

The 4-D lattice approach meets two problems which 
made the 31? approach so attractive. At first, only the 
bosonic sector can be studied because of the problems to 
deal with chiral fermions. Second, simulations are difficult 
due to the presence of different length scales. The inverse 
Matsubara frequency, 1 / (2ttT) , is much shorter than the 
two correlation lengths, 1/mjy (T) and l/mjy(T), with the 
screening masses differing strongly from the zero temper- 
ature masses Mw and Mh- If the temporal lattice extent 



N T should be not too small, they can be accommodated on 
a 4D lattice only if this is anisotropic. This requires more 
couplings to be tuned than in the isotropic case when one 
first has to match the Mh/Mw ratio and the renormalised 
coupling 54(^4) on a T — (N T — N s ) lattice. 

The 3D approach represents itself as a radical solution 
to both problems. It is based on perturbative dimensional 
reduction in the continuum combined with lattice pertur- 
bation theory ||, leaving to lattice simulation only the 
non-perturbative description of the softest modes whose 
behaviour governs the transition. All non-zero Matsub- 
ara modes are integrated out leaving an effective, super- 
renormalisable theory in 3-D with the smallest length scale 
l/(27rT) removed. Chiral fermions can be implicitly dealt 
with at this step. In a second step the next heaviest modes 
related to the Debye mass mjj ~ g T (the Aq components 
of the gauge field) are integrated out. Thus an effective 
3-D S't/(2)-Higgs theory emerges with the action 



S3 



(1 



A, 



(1) 



It has dimensionful, renormalisation group invariant cou- 
plings g| and A3 and a running mass squared m^/^). It is 
put into correspondence to a lattice model with the action 
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The lattice couplings are (with a suitable parameter M^) 
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which can be expressed in terms of 4D couplings and 
masses. The parameter is approximately equal to the 
zero temperature physical Higgs mass. The summation in 
(H) is taken over plaquettes p, sites x and links I — {x, /i}. 
The gauge fields are represented by unitary 2x2 link ma- 
trices U~x.fi and Up denotes the SU(2) plaquette matrix. 
The Higgs field is parametrised as follows: <P X = p x V x , 
where p\ = ^tr(<P^<P x ) is the Higgs modulus squared, and 
V x is an element of the group SU(2). 

The bare mass squared is related to the renormalised 
m§(/i3) (we choose /13 =33) through a lattice two-loop 
calculation ]l(| giving m|(^ 3 ) 



l l-loop 



l 2-loop- 



The lattice results obtained in the 3D approach indicate 
the validity of the dimensional reduction (without addi- 
tional operators in the action) near to the transition tem- 
perature for Higgs masses in the range between 30 and 
240 GeV. The relations between the parameters of the ef- 
fective 3D model and the physical quantities Higgs mass 
and temperature are derived in ||. 

According to recent lattice studies |4|,||,|ll} the stan- 
dard electroweak theory ceases to possess a first order 
transition for a Higgs mass Mh > 72 GeV. Therefore, 
taking the newest lower bound of the Higgs mass into ac- 
count Jl^l , the standard model does not pass through a 
true phase transition at the electroweak vector boson mass 
scale. This and the small amount of CP violation in the 
standard model seem to rule out the possibility to explain 
the BAU generation without new physics. Therefore the 
phcnomcnological interest has moved to extensions of the 
standard model, with minimal supersymmetric extensions 
(MSSM) being the most promising variant. First interest- 
ing lattice results have been published recently [ p^[ . 

From the point of view of non-perturbative physics in 
general, the lattice version of the standard Higgs model 
is still interesting as a laboratory for investigating the be- 
haviour of hot gauge fields coupled to scalar matter, for 
the characterisation of possible bound states, for the un- 
derstanding of real time topological transitions and the 
role of embedded topological defects (m) at the transi- 
tion. It has been important as a cross-check for analytical 
approximation schemes and will be so in future. 

We are concentrating in this paper on the qualitative 
change of the spectrum of screening states that happens 
across the phase transition and what remains of it in the 
crossover region (at somewhat higher Higgs mass). This 



would not be possible without a systematic evaluation of 
the wave functions in configuration space that, for the first 
time in this context, is attempted here. Some intermediate 
results have been published before JL5| ]. 

The rest of this paper is organised as follows. In Sect. || 
we discuss the cross correlation technique and the operator 
set as we have used them. The wave functions and masses 
of the ground state and excited states near to the endpoint 
of the phase transition are presented in Sects. 3 and 4 for 
the symmetric and the Higgs phase at = 70 GeV, re- 
spectively. This complements a recent study in fl6| which 
has been performed, however, at much lighter (M^ = 35 
GeV) Higgs mass. Section 5 contains our results for the 
dynamics of the spectrum with increasing temperature in 
the crossover region for Mjj = 100 GeV, slightly above 
the endpoint of the first order transition, which can be 
compared with results of Ref. JlTj obtained at a markedly 
larger Higgs mass Mjj « 120 GeV. In Sect. 6 we sum- 
marise our results. The Appendices contain tables of the 
measured mass values obtained from this analysis and de- 
tails how to construct the operator set for the quantum 
numbers in 2 + 1 dimensions on the lattice. 



2 The Cross Correlation Technique in the 3D 
Model 

To study simultaneously the ground state and excited 
states (as well as their wave functions) one has to con- 
sider cross correlations between (time-slice sums of) oper- 
ators Oi from a complete set in a given J PC channel with 
quantum numbers J (angular momentum) , P (parity) 
and C (charge conjugation). According to the transfer ma- 
trix formalism, one should be able to write the connected 
correlation matrix at time[j] separation t in the spectral 
decomposition form 



with 



<vac|Oi|* (n) > 



(4) 



(5) 



where | is the n-th (zero momentum) energy eigen- 

state. The vacuum state is dropped from this sum due 
to the connectedness of the correlator. By suitable diago- 
nalisation this allows to find masses and wave functions of 
the lowest mass screening states (ground state) and higher 
mass excited states in the various J PC channels. 

However, in practice one has to choose a truncated 
set of operators Oi, (i — 1, . . . , N). The hope is that the 
lowest lying states (k — l,...,r, with r << N) might 
not be essentially affected by the truncation and that 



in the 3D approach "time" is x$ 
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their masses and wave functions can be approximately ex- 
tracted from the direct eigenvalue problem for C'ij(t) 

^^(*)-^ n) =A(")(t)^. (6) 
j 

Experience shows that this solution suffers from big sys- 
tematic truncation errors. Solving instead the generalised 
eigenvalue problem 

£C,y(f) *f 5 = A(")(Mo)£c^ ) fj n) (7) 
j j 

or 

^■(Mo^j-^A^o)^, (8) 

3 

with 

Cy(Mo) - (c-i(i )C(t)C7-'(to)) , (9) 

(t > to, where practically t — 0,1,2,..) errors related 
to this truncation can be kept minimal |t8|[nj. Practi- 
cally the decomposition of the matrix C(tt ijis performed 
using a Cholesky decomposition C(to) = LL T . The re- 
maining problem is that of diagonalizing a symmetric ma- 
trix C(t,ta) = L~ 1 C(t)L T _1 with rotated eigenvectors 

The optimised eigenfunctions \p( n > in the chosen op- 
erator basis (obtained with a small distance to) give an 
information about the overlap of the source operators Oi 
with the actual eigenstates \&(™'). Due to truncation the 
eigenvectors of (Q) are not supposed to be orthogonal to 
each other since C~ 1 (to)C(t) is not a symmetric matrix. 

The components of the eigenvectors K™' (n < N) are ef- 
fected by terms of the order 0(exp(— mN+it)), where TV is 
the number of used operators. Only in the limit of a com- 
plete set of operators the eigenvectors become orthogonal. 
Therefore, the (non) orthogonality between different states 
provides a criterion for the completeness of the operator 
set. 

The masses rrv- n > of these states are obtained by fitting 
not the eigenvalues A^™^ of (||) but the diagonal elements 
//("), 

^(t,t )=J2 <3tf(Mo)#j n) , (10) 

to a hyperbolic cosine form with t in some plateau region 
of a local effective mass that is defined as 

my t, to) = log , . , -. (11) 

off v ' ' 6 //(™)(t + l,t ) v y 

The wave function components with respect to the op- 
erator basis characterise the coupling of source operators 
Oi to the bound states (lowest mass or excited) in the J 
channel. Using operators of different extension transverse 



to the correlation direction a spatial resolution of the op- 
timised wave function can be achieved. Practically only a 
subset of operators with fixed quantum numbers can be 
included. The admitted operator set includes gauge invari- 
ant operators properly chosen with respect to the lattice 
symmetry and quantum numbers. For completeness, the 
choice of the gauge invariant operators in 2 + 1 dimensions 
on the cubic lattice corresponding to the quantum num- 
bers mentioned is discussed in Appendix |A|. Though the 
results might be in principle known a reasonable presen- 
tation for the lattice S , J7(2)~Higgs model similar to that 
of Ref. |23| for 3+1 dimensions was not available. 

In order to associate a spatial structure to the states 
under study one has to use operators which correspond 
to various extension. One way is to define a set of opera- 
tors using smearing techniques in several variants. In the 
present context this has been practised in Refs. |p0j|l6l , 
|l7|| , smearing gauge links and Higgs fields. The smearing 
parameter determines a well-defined mixture of various 
lengths in the source operator and has to be optimised by 
a variational procedure (strictly speaking, for all states 
separately). 

In contrast to this smearing technique, we have chosen 
the other extreme and collected only a few types of opera- 
tors Oi in our base (properly chosen with respect to lattice 
symmetry and quantum numbers) but with a wide span 
of sizes I in lattice spacings. Such a basis allows to obtain 
information on the spatial extension of a bound state with- 
out going through a variational procedure (for the lowest 
and the excited states). In principle, smearing the link 
and Higgs variables underlying these operators would be 
conceivable for further optimising the states under study. 
Here, we have restricted ourselves to Higgs strings and 
Wilson loops with varying spatial extent. 

In the three quantum number channels investigated 
here we have admitted the following operators (with /i = 3 
reserved for the correlation direction, see Appendix |a|): 

++ : Pi 

S xA (l) + S X>2 (1) 

W x ,i,i(l) + W x , 2 ,i(l) (12) 
(quadratic Wilson loops of size I x I) 

1~ : V x b A (l) + V x b t2 (l) 

2 ++ : S xA (l) - S x , 2 (l) 

Here the following notation is used (I denotes the string 
length in lattice units) 

SxA 1 ) = l^tUx-n ■ ■ ■ U x+{ i-i)^ x+ ip) , (13) 
V x b Jl) = hr(r b <PtU X4l . . . U BS+{l _ 1) M&*+ip) • (14) 

Since already the masses of the lowest states in the 2 ++ 
channel are relatively heavy we have not included Wilson 
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loop contributions here. Therefore, the wave functions of 
2 ++ excited states have to be taken with care. 

~ (n) 

In our procedure, finding the eigenfunctions W- is 
tantamount to determine an optimal source operator 
for the eigenstate |\I>( n )) as a superposition of the original 
operators Oi, 

N 

0(»)=^^")O i! (15) 

where the coefficients iP^ are related to the normalised 

solutions <Pj of the generalised eigenvalue equation (||) 

by <?f l) = Lj. -^W/diWfe.io)) 1 ' 2 . These coefficients 
define the weight of the chosen operators contributing to 
G>(") and express the spatial extendedness of the eigen- 
state |<3>(™'} under study. 

In the 1" and 2 ++ channels the index i is used to la- 
bel the size I of the corresponding Higgs string operator in 
lattice units. In the case of the ++ channel the label I = 
refers to the operator p^,, and for I > 1 we have to distin- 
guish between the I th contribution from Higgs strings and 
that from Wilson loops. 

We have performed simulations using the update algo- 
rithms as described in our previous works J3| . In the error 
analysis of the wavefunctions and masses we have used the 
jackknife technique. As an example of the outcome of the 
analysis we show in Fig. Ill the exponential decay of the 




5 10 15 

t 



Fig. 1. Example of the exponential decay of p} n \t,to) for the 
five lowest states 



diagonal elements of the correlation matrix /i( n \t,to) for 
the first five states (n = 1, . . . , 5) in the ++ channel as 
they are found in the symmetric phase (at (3q — 12, 0h = 
0.3434) for a Higgs mass of M* H = 70 GcV. The mass 
estimates are obtained by fitting p^ n '(t,ta), throughout a 
plateau region of the effective mass, to a hyperbolic co- 
sine form, taking into account the estimated errors of the 



0.20 -° 
0.18 - 



a o.i6- 
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Fig. 2. ++ ground state effective mass obtained from cross 
correlations (o), compared to results solely from Sx,i+S x ,2 (□) 
and pi (-*-) . 



correlation functions. The improvement achieved by using 
cross correlations is demonstrated in Fig. |[ As expected, 
the cross correlation method leads to lower mass values. 

The measurements for a Higgs mass = 70 GeV 
(A3/33 ~ 0.09570) were performed below but near to the 
endpoint of the phase transition. At both sides, in the sym- 
metric and in the Higgs phase, we used a respective line 
of constant physics corresponding to m\{g1)/g\ w 0.00023 
in the symmetric and rri^(g^)/gl w —0.043 in the Higgs 
phase. This corresponds to T « 1.02T C and T w 0.98T C , 
respectively [|[||] . The gauge coupling ftc has been varied 
from 8 to 16 in order to examine the extrapolation to the 
continuum. The choice of the actual simulation parame- 
ters as near as possible to the phase transition line was 
dictated by the necessity to avoid tunnelling between the 
phases. Besides the used lattice size of 30 3 we have addi- 
tionally investigated the wave functions in the symmetric 
phase for the 1 and 2 ++ channels on a 50 3 lattice which 
permitted to increase the range of I in the operator basis. 
On a lattice of given size, the spatial extension of opera- 
tors is restricted to half of the lattice size. 

At the larger Higgs mass of M* H = 100 GeV (A3/5I rj 
0.1953), above the endpoint of the phase transition (J, 
||, we have studied the spectral change with decreasing 
temperature (increasing (3h), from the "symmetric" to the 
"Higgs" side of the crossover line, only at one fixed (3q = 
12. 

The accumulated statistics of our measurements in the 
symmetric and Higgs phases as well as in the crossover re- 
gion for each pair of /3q and [3h values is summarised in 
Table [lj By measuring the autocorrelations for the oper- 
ators of interest we determined the frequency of measure- 
ments during the updating in order to retain only inde- 
pendent configurations. 
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Table 1. Used statistics per (/3g, Ph) in the measurements 



Phase 


Channel 


GeV 


L 3 


Pa 


independent 
configurations 


symmetric 


1 — ,2++ 


70 


50 3 


8,12,16 


20000 




0++.1— ,2++ 


70 


30 3 


8,12,16 


12000 


Higgs 


0++1— ,2++ 


70 


30 3 


8,12,16 


12000 


crossover 


++ ,l~ ,2++ 


100 


30 3 


12 


4000 















3 Bound States in the Symmetric Phase at 



M 



H 



70 GeV 



Using the cross correlation technique we were able to ob- 
tain the wave function squared corresponding to the opti- 
mised operator for each individual state in the spectrum. 
Being functions of a physical distance, the squared wave 
functions are shown immediately vs. I a g\ in order to 
overlay data from measurements at various gauge cou- 
plings (lattice spacings) taken along a line of constant 
physics. 

The results for the ++ channel are collected in Figs. Il- 
ls] for the squared wave functions. The contributions from 



1.2 


M|,l, 


IIM I 1 1 1 1 1 1 [ 


1.2 


.'1 1 1 


1 1 1 1 1 '_ 


1.0 




□ /3 G = 8 

O &=12 : 


1.0 




□ fc=8 






a 3,-ie 




" 1 


a G =16 


o.a 






0.8 


o 












- k 












~ 




0.6 






0.6 














1 




0.4 






0.4 


: '§ 












"8 
















0.2 






0.2 


% 












w 











I o 


.1 1 1 


^%*>0O*>^ □ y □ y a- 







2 4 6 




2 


4 6 






1 o gf 






I ° 9l 



Fig. 3. Squared wave function of the ground state in the ++ 
channel, measured on a 30 3 lattice in the symmetric phase; 
left: W B ,i, 2 (i) + W.,a,i(0. ri S ht: S..i(0 + S -M; 1 = 0,..., 15 



the Higgs string and Wilson loop operators are shown sep- 
arately in order to identify clearly Higgs and H^-ball exci- 
tations. We observe no mixing of these two operator types 
in the Higgs ground state and the first excitation. The 
second excited state in this channel consists of a pure ex- 
citation of gauge degrees of freedom (d.o.f.) and can be 
identified with a VF-ball in analogy with the glueballs of 
pure SU(2). Our results on this decoupling confirm the 
observations in jl6| made at a much lighter Higgs mass in 
the symmetric phase near to the strongly first order phase 
transition. 




Fig. 4. Same as Fig. H for the first excited state 



S G =1 2 
? G =16 



if H 



O (3 E =12 



Fig. 5. Same as Fig. ^ for the second excited state 



The squared wave functions for ground and excited 
states in the 1 and 2 ++ channels are presented in Fig. || 
and @. 

At this point it is due to come back to the issues of 
completeness (of the operator set) and orthogonality (of 

the states tyj™^). We have already noticed that the eigen- 
vectors of (pj) become orthogonal only in the limit of a 
complete operator set, and that the scalar product be- 
tween different states should be used to examine the de- 
gree of completeness. This test shows what are the dif- 
ficulties in practice. We have studied the scalar product 
between the ground state and the first excitation as a func- 
tion of the maximal length l max of the operator set S X ^{1) 
by gradually clipping the set of operators from / ma ic=25 
(on the 50 3 lattice) down to zero. One would expect that 
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this function tends to zero for L 



co. Actually the 



scalar product reaches a minimum at some finite l m ax and 
starts to grow at bigger maximal length due to statistical 
and numerical errors in the case of larger operators and 
matrices. 

Inspecting the wave function of the ground state in 
the 1~~ channel we observe that for very large operators 
the contribution to this state vanishes. This is true to a 
high accuracy for the largest physical volume that we have 
considered, at Pa — 8. If we chose higher (3g (smaller 
volumes) along the line of constant physics (in order to 
explore the approach to the continuum) the operator set 
in use with fixed l ma x = L/2 does not describe anymore 
the whole wave function. There are fake contributions ac- 
cumulating in the contribution of the longest operator in 
the set which, in principle, belong to more extended oper- 
ators which are not included. A similar behaviour can be 
observed for the excited states if one compares the wave 
functions varying the operator set. For the second excited 
state shown in Fig. ^ it is clearly seen that in the case 
of /?g=16 (small physical volume) this state is not satis- 
factorily described by the limited operator basis. An opti- 
mal continuum limit for the wave function would require 
to choose the same maximal operator length in physical 
units, i.e. doubling the length in lattice units when one 
goes over from Pq = 8 to Pq = 16. This is, however, 
difficult to accomplish with restricted resources, not only 
because of the lattice size. Larger cross correlation ma- 
trices would additionally require larger statistics to get 
stable results from the diagonalisation procedure. 

Considering only wave function data obtained at fixed 
Pg (tolerating a non-perfect overlay of data from differ- 
ent a in a physical scale) the presented figures illustrate 
perfectly that the number of zeroes of the wave functions 
corresponds to the level of excitation (knot rule). 

Due to the mentioned difficulties, it becomes more and 
more difficult to extract the masses of higher excitations, 
in particular if already the ground state is very heavy and 
if the lattice becomes dangerously small (when we attempt 
to study the continuum limit) compared to the physical 
extension of the wave function corresponding to the state 
under discussion. 

For the ground states we are now able to reexamine our 
previous spectrum investigations [|j where we measured 
the mass without use of the cross correlation technique. 
In the ++ channel the agreement is satisfactory since 
the ground state wave function is well dominated by the 
shortest operators like {l/2)tv{^U x ^^ x+ p). Mass mea- 
surements in the 2 ++ channel and for the VK-ball states 
have proved to be difficult in our previous studies, and 
we did not report on that in Ref. |J. From the present 
analysis we recognise that this difficulty is related to the 
large spatial extendedness of these states. 

Using the cross correlation technique we have found 
in the symmetric phase the ground state mass and the 
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Fig. 7. Squared wave function in the 2 channel, measured 
on a 50 3 lattice in the symmetric phase; ground state and first 
excited state 



first excitations. We restricted ourselves to an identifica- 
tion only of those masses with a still reasonable plateau 
behaviour and a not too large mass value in lattice units. 
Notice that the operators are classified in the angular mo- 
mentum J only modulo 4; therefore we cannot exclude 
that states with higher masses have got contributions from 
states with higher J. 

The masses of the first excitations in the symmetric 
phase (for 2 ++ only of the lowest state) are presented in 
Fig. H and the fitted masses are collected in Tables |9- 10 of 
Appendix [b|. We tried to illustrate whether an approach 
to the continuum limit is already indicated by the present 
measurements. The big errors of the excited state in the 
1 channel and, even more, of the ground state in the 
2 ++ channel show where the problems are. 

Concluding this Section we would like to stress that 
qualitatively there is no difference in the spectrum on the 
high temperature side of the phase transition for light or 
relatively heavy Higgs masses as long as the phase tran- 
sition persists. The gauge d.o.f. decouple from the Higgs 
excitations and form VF-balls similar to confined 31? pure 
SU (2) gauge theory M . 



4 The Spectrum in the Higgs Phase near to 
the Endpoint of the Transition 



On the Higgs side of the phase transition pure gauge mat- 
ter (W-ball) excitations are not expected to be present in 
the spectrum. This hypothesis has been numerically ver- 
ified at light Higgs masses @. Now, at = 70, we 
observe in the Higgs phase near to the transition the fol- 
lowing. In the 1" and 2 ++ channels the spectrum looks 
similar to that in the symmetric phase. This was expected 
from earlier studies for the 1 channel. Therefore the 
squared wave functions are not shown explicitly. 
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Fig. 6. Squared wave function in the 1 channel, measured on a 50 3 lattice in the symmetric phase; ground 
state and two excited states 
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Fig. 8. Mass spectrum in the symmetric phase as function of 1//3g 



Note, that in the 2 ++ channel we have not used oper- 
ators with pure gauge degrees of freedom. Therefore, we 
cannot observe the expected difference for 2 ++ VF-balls.^] 

In the ++ channel, however, our operator set is suf- 
ficient to observe a marked difference between the phases 
which is not in accordance to naive expectations. As a 
characteristic feature we observe the mixing between the 
two operator types, W-ball operators (pure gauge d.o.f.) 
and operators projecting onto Higgs states. The squared 
ground state wave function looks similar to that in the 
symmetric phase, contributions from Wilson loops to the 
operator projecting onto the lowest mass states are absent 
(Fig. ||). However, already the first excited Higgs state con- 
tains a noticeable contribution from Wilson loop operators 
of almost 20 percent (Fig. |To|) . Earlier measurements at a 
Higgs mass of — 35 GeV Jl7j (where the phase tran- 
sition is very strong) did not indicate such a mixing. We 
interpret this mixing of Higgs and gauge d.o.f. as a signal 
of the near endpoint of the phase transition. Deeper in 
the Higgs phase the contribution from gauge degrees of 
freedom is expected to disappear also at this high Higgs 



O 00=12 
« fo-16 




Fig. 9. Squared wave function of the ground state in the ++ 
channel, measured on a 30 3 lattice in the Higgs phase; left: 

W«,i,2(I) + W«,3,i(i). ri s ht: &,i(0 + SxAi); i = o, • • • , 15 



Results for the W-balls with J — 2 can be found in 



mass. This tendency has been checked in our simulations 
at Mjr = 100 GeV to be discussed in the next section. 
This suggests a scenario according to which the first and 
second excited states present on the symmetric side (Higgs 
and H'-ball, respectively) merge into one common state on 
the Higgs side. The remaining gauge degrees of freedom 
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Fig. 10. Same as Fig. g for the first excited state 

are fading away from that state with further lowering of 
temperature (increasing Ph)- 

The second excited state in the ++ channel on the 
Higgs side (Fig. [ll]) is qualitatively a state with Higgs 




Fig. 11. Same as Fig. ^for the second excited state 

excitations (no VF-ball contributions) followed by a third 
excited state which looks dominantly like a VF-ball.[] 

In Fig. [l2| our results for the masses at three differ- 
ent gauge couplings on a line of constant physics are col- 
lected, the corresponding values can be found in Tables 
|Tl| and Similar as before, the lowest mass state has a 
reasonable continuum limit whereas with increasing mass 
the accuracy rapidly deteriorates and an extrapolation to 
vanishing lattice spacing becomes difficult. 



5 The Spectral Change at the Rapid 



Crossover for M 



H 



100 GeV 



We have noticed that the spectra of the two phases dif- 
fer mainly with respect to the contributions of the Wil- 
son loops (gauge d.o.f.) to the excited states in differ- 
ent channels. Therefore we study this change in more de- 
tail while continuously passing the crossover line (chang- 
ing the hopping parameter Ph) at fixed gauge coupling 

3 Its wave function is not shown here. 



(3q — 12. Sufficiently above the phase transition endpoint 
(for large enough scalar self-couplings) no large autocor- 
relation times are expected which would prevent us to 
cross that region of still rapidly changing thermodynami- 
cal observables and to determine the physical excitations. 
As long as a true phase transition exists, due to tunnelling 
one would measure in the ++ channel not the actual low- 
est mass in the spectrum related to the respective phases 
but rather a characteristic correlation length characteris- 
ing the transition itself. 

Our simulation have been performed at a scalar-gauge 
coupling ratio X^/gl « 0.1953 corresponding to the ap- 
proximate Higgs mass parameter M H = 100 GeV. We 
have used the same operator set as in the previous analy- 
sis. The statistics of 4000 independent configurations per 
Ph value limits the capability to analyse higher excita- 
tions. Nevertheless, we found a behaviour very similar to 
our results obtained at M H = 70 GeV. 

The similarities concern both the high temperature 
side of the crossover (where one expects thermodynamic 
properties being close to those of the symmetric phase 
at smaller Higgs mass) and the region very near to the 
crossover line on the so-called Higgs side of the crossover 
(resemb ling the Higgs phase at slightly lower Higgs mass) . 
In Figs. yj| we present the spectrum of the lowest states in 
the ++ , 1" and 2 ++ channels as function of m^gf )/g| 
(or Ph) over a wide interval above and below the crossover 
(i. e. in temperature). 

By inspecting these Figures we localize the crossover 
at wKffD/gf ~ —0.05. Of particular interest for us is 
to study the mixing of Higgs and gauge degrees of free- 
dom in the ++ channel. Looking at the excited states 
in this channel of Figs. [l3] we conclude that the scalar 
and gauge sector are approximately decoupled as long as 
one keeps away from the crossover line on the high tem- 
perature side. This has already emphasised in Ref. [[l6| 
[l7fl where cases of smaller (characterized by the existence 
of a strong first order phase transition) and larger Higgs 
self-couplings (far above the endpoint of the phase tran- 
sition) have been studied. The mass of the lowest IF-ball 
state (full triangle) is roughly independent of Ph as long 
as one does not come too close to the crossover. Thus, on 
the high temperature side of the crossover, the ordering 
of states qualitatively resembles the spectrum at smaller 
values of Higgs self-coupling (at M H = 70). If one ap- 
proaches the crossover temperature the mass of the first 
(Higgs-like) excitation is moving up towards the lowest 
VF-ball state whose mass decreases. At some point we 
observe a growing admixture to the Higgs excitation by 
contributions from Wilson loop operators which reaches 
about 17 %. Then the wave function of the remaining 
state looks similar to figure [h]. At higher Ph (lower tem- 
perature, deeper in the would-be Higgs phase) the ad- 
mixture from pure gauge d.o.f. disappears again from this 
state. This has been explicitly checked for the first ex- 
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Fig. 12. Mass spectrum in the Higgs phase as function of 1//3g 



Fig. 13. Spectra near the crossover in the channels ++ , 1 and 2 ++ ; in the ++ channel open triangles denote 
Higgs states, full triangles W-ball states and open circles Higgs states with an admixture of excited gauge d.o.f. 



cited state in the ++ channel corresponding to the data 
point at /3h = 0.37, fn\{g1)/g^ ~ —3.4 in table [l5| (far on 
the Higgs-side) where we did not observe any contribution 
from Wilson loops. 

6 Summary and Conclusions 

In this paper we complemented our numerical study of 
the electroweak phase transition using the 3D SU(2)- 
Higgs model by an investigation of the ground state and 
some excited states in three J PC channels. Our interest 
was focused on Higgs screening masses slightly below and 
above the end of the first order phase transition where we 
wanted to understand the qualitative relations between 
the higher and lower temperature regions. This investi- 
gation was done with the help of cross correlation matrix 
functions between different operators with the same quan- 
tum numbers. We were able significantly to improve the 
old ground state mass estimates which were obtained just 
by fitting a single correlation function over the plateau re- 
gion of the effective mass. In our earlier work, the proper 
extension of the operators had been selected by hand, 



according to a signal vs. noise criterion for various di- 
agonal correlators. Implicitly this procedure, too, gives 
some information concerning the best-projecting operator 
that relates the vacuum state to the lowest mass bound 
state in the given channel. Today, many techniques are 
known allowing to construct source operators with im- 
proved projection properties. For the SU (2)-Higgs and 
the pure SU (2) model these are smearing techniques which 
were investigated and used in our present context in Rcfs. 
[^0|[l(|[l7]]. These improved operators take the extended 
character of the states into account, too, although in a 
summary way. 

Using the cross correlation technique the extended char- 
acter can be specified in a more detailed and systematic 
way. Within a preselected set of operators, the method 
finds the projection strength of a each operator onto the 
first few states with quantum numbers compatible to the 
operators admitted. The result can be represented as a 
wavefunction. Moreover, different Fock space components 
(Higgs bound states and gauge-ball states) can be treated 
in parallel. In the present paper this feature has been 
proven essential for the understanding of excited states 
and to point out similarities and differences between neigh- 
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bouring regions of the extended phase diagram. Our oper- 
ators were classified in a very simple way (extended Higgs 
strings and Wilson loops without internal structure) with 
a clear reference to the corresponding length in lattice 
units. This would allow us to interpret the results in terms 
of a configuration space wave function. Anticipating this 
possibility and just to make the first step, our operator 
choice was the simplest one to incorporate the notion of 
spatial extension. It was possible to see the particular ef- 
fects that generically occur when the admitted operator 
set becomes too small (in attempts to go too far towards 
the continuum limit a — > 0) to cover the state under dis- 
cussion. The continuum limit of the screening masses and 
wave functions has to be accompanied by correspondingly 
larger lattices with the same physical volume. But then, 
simply increasing the operator set simultaneously, would 
require much higher statistics to reach the necessary accu- 
racy and stability in the diagonalisation procedure. This 
is a reason to include, in a next step of improvement of 
the method, smearing of the fields entering the operators 
and/or the construction of blocked operators. 

Our first measurements presented here were done at a 
Higgs mass of = 70 GeV, near to but still below the 
endpoint of the transition. We found that the physics at 
this point is already influenced by the endpoint. The W- 
ball states characteristic for the symmetric phase do not 
disappear in the Higgs phase as long as one keeps near to 
the transition line. However, the first VK-ball state appears 
here as an admixture to the first excited Higgs state. A 
second PT-ball state could also be seen in the Higgs phase 
very near to the phase transition. 

Our measurements in the crossover region have shown 
that also there the VF-ball state will finally disappear 
deeper on the low temperature side of the crossover. In our 
simulations at = lOOGeV we could study the mixing 
of the first W-ball with the first Higgs excited state, both 
being well separated on the symmetric side of the crossover 
(similar to the symmetric phase at lower Higgs mass), in 
some detail by continuously monitoring the spectral evolu- 
tion across the crossover line. Over most of the parameter 
space on the high temperature side of the crossover re- 
gion, however, decoupling of the gauge and scalar sector 
has been established. 



A Angular Momentum in 2+1 Dimensions 

Gauge invariant operators can be constructed in the 3D 
SU (2)-Higgs model and can be used to define states and 
wave functions. They are classified with respect to the 
quantum numbers angular momentum J, parity P and 
charge conjugation C. While the classification according 
to C and P is trivial (the eigenvalues of C and P have 
the same sign for operators which live in the same x — y- 
planc in 2 + 1 dimensions) the classification according to J 



requires some care due to the discreteness of the rotation 
group. The results are known, we collect here the main 
results for the convenience of the reader and to apply them 
to the operators of interest. For the 3 + 1 dimensional case 
the classification can be found in Refs. Some of 

the results are already described in pffi we follow the 
description in f25|| . 

To define the angular momentum of a given operator 
on the lattice the transformation properties with respect 
to the discrete lattice symmetries have to be studied. In 
our case the continuous symmetry group is SO (2) which 
gets restored in the continuum limit. In that case the an- 
gular momentum is represented by the irreducible repre- 
sentations (IR) of SO(2). On the cubic lattice the group of 
two-dimensional rotations and reflections is the dihedral 
group P4. This group is non-Abelian and has 8 elements 
(d = 8) and 5 irreducible representations. The conjuga- 
tion classes C with their number of elements nc and the 
symmetry transformations are collected in Table 0. The 



Table 2. Symmetry transformations of T>4 group 



C{n c ) 


angle 


1(1) 


Identity 


cf(i) 


rotation (0,0) <j> — tv 


C±{2) 


rotation (0,0) <j> = ty/2 


C a (2) 


reflection x, y 


C£(2) 


reflection diagonal 



corresponding IR's of V4 are denoted by A\, A2, B\, E>2 
and E, the first four are one-dimensional, the last one has 
dimension 2. 

The characters of P4 are computed by reducing the 
characters of the full SO (2) group (including reflections) 
(table H) to those of the subduced representation Dq by 
taking the discrete angles of table ^. The result is shown in 
table 0. The characters of the conjugation classes of V4 in 

Table 3. Characters x J of full SO (2) symmetry (C(4>)) in- 
cluding reflections (cr„) for various angular momenta J 








1 


2 


3 


4 


1 


1 


2 


2 


2 


2 


cy>) 


1 


2cos<^> 


2cos2c/> 


2cos30 


2cos4<^> 


a v 


1 















the IR of this group are taken from the character table of 
[p5| and are shown in table ||. Using the characters given 
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Table 4. Characters \c of conjugation classes C in the sub- 
duced representation Dq 



C\J 





1 


2 


3 


4 


1(1) 


1 


2 


2 


2 


2 


Cl(l) 


1 


-2 


2 


-2 


2 


C 4 (2) 


1 





-2 





2 


C a (2) 


1 














a (2) 


1 















Table 7. Characters and multiplicities for the operators 

s.M and v£, M (J) 



xb R 


i(i) 


cl(i) 


C 4 (2) 


Ca(2) 


C 2 (2) 


Sx,fj, (0 


2 


2 





2 







2 


-2 













Ai 


A 2 


Si 


B 2 






1 





1 






















1 



Table 5. Characters \c °f the conjugation class C in the IR 
of ©4 



C\IR 


Ai 


A 2 


Bi 


B 2 


E 


1(1) 


1 


1 


1 


1 


2 


Cf(l) 


1 


1 


1 


1 


-2 


C 4 (2) 


1 


1 


-1 


-1 





C 2 (2) 


1 


-1 


1 


-1 





C 2 (2) 


1 


-1 


-1 


1 






Table 6. Multiplicities rrij of IR of T>4 in the subduced rep- 
resentation Dq for angular momentum J 



J\IR 


Ai 


A 2 


Bi 


B 2 


E 





1 














1 














1 


2 








1 


1 





3 














1 


4 


1 


1 












in tables ||| and || we can compute the multiplicities for the 
IR of the dihedral group T>± for angular momentum J: 



JR 



ncXc R Xc » 



(16) 



c 



they are shown in table ||. 

The multiplicities for a given operator are obtained 
by inspecting its properties under those discrete transfor- 
mations. The extended operators (given in ji~3|)) trans- 
form as vectors with positive parity (S Xjfl (l)) or with neg- 
ative parity (V^Al)). The characters and the multiplici- 
ties are shown in table ?]. We observe that the operator 
S x ,n(l) transforms simultaneously under the IR Ai and 
B\, whereas V x „(?) transforms with E only. 



Using a projection operator P IR one is able to project 
onto one irreducible representation in the case of S x ^(l). 
A projection operator for a given IR is constructed as 
the sum over all conjugation classes C weighted by the 
character Xc* for that IR (table ||) : 



P 



IR 



(17) 



c 



The C have to be taken in a matrix representation taking 
into account parity. 

For the operator S XtfM (l) the projection matrices can 
be used in the form: 



P- 



1 1 
1 1 



P 



1 -1 
-1 1 



(18) 



From here the fixed angular momenta of the operators 
to be used are defined as shown in table @. It is obvious 



Table 8. Quantum numbers and gauge invariant lattice oper- 
ators 



0++: 


S x ,i(l) + S Xj2 (l) 




W X ,l, 2 (l)+Wy,l,2(i) 


1~: 


<!(!) + Vtf, 2 (l) 


2++: 


S x ,i{l) - S x , 2 (l) 




W^,i,a(i)-Wi,,i,a(i) 



from table ^ that all operators describing a J = angu- 
lar momentum state also overlap to a J = 4 state. These 
higher states are expected to be sufficiently heavier than 
the lightest J = state. Similarly this also holds for the 
other operators. Thus for the cubic lattice only operators 
with J(mod 4) can be constructed what restricts the iden- 
tification of highly excited states in the channels discussed 
here. 
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B Tables of Masses for Ground States and Excitations 



Table 9. Masses of ++ states in units of gl at M* H = 70 GeV in the 
symmetric phase (30 3 lattice) 



Pg 


Ph 


0++ 


0++* 


0++ 

U W-ball 


Q + +** 


0++* 

U W-ball 


8 


0.3490 


0.504(16) 


1.628(64) 


1.998(37) 






12 


0.3434 


0.398(18) 


1.515(33) 


1.902(26) 


2.679(75) 


2.868(123) 


16 


0.3407 


0.404(59) 


1.824(100) 


2.080(83) 


3.604(172) 


2.736(228) 



Table 10. Masses of 1 and 2 ++ states in units of gl at M H = 70 GeV 
in the symmetric phase (30 lattice) 



Pg 


Ph 


1"" 


1—* 


2++ 


8 


0.3490 


1.122(35) 


2.054(87) 


1.498(68) 


12 


0.3434 


1.074(14) 


1.857(68) 


1.518(78) 


16 


0.3407 


0.992(68) 


2.028(259) 


1.788(21) 



Table 11. Masses of ++ states in units of gl at M* H = 70 GeV in the 
Higgs phase (30 3 lattice) 



Pg 


Ph 


0++ 


++ * + gauge d.o.f. 


++ 

u W-ball 


Q+ + ** 


8 


0.34970 


0.380(4) 


1.836(36) 


2.296(148) 


3.206(519) 


12 


0.34368 


0.390(3) 


1.815(90) 


2.862(71) 


3.009(126) 


16 


0.34088 


0.404(4) 


2.268(71) 


3.308(244) 


3.668(112) 



Table 12. Masses of 1 and 2 ++ states in units of gl at M H = 70 GeV 
in the Higgs phase (30 3 lattice) 



Pg 


Ph 


1"" 


l — , 


2++ 


2++* 


8 


0.34970 


0.820(8) 


1.894(134) 


1.764(184) 


2.422(288) 


12 


0.34368 


0.798(6) 


1.767(189) 


1.830(48) 


2.811(63) 


16 


0.34088 


0.820(12) 


2.420(151) 


1.940(100) 


2.944(344) 
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Table 13. Masses of ++ states in units of g\ at M* H = 100 GeV 
(crossover region) depending on (5h 





Ph 


0++ 


0++* 


0++ 

U W — ball 


symmetric 


0.34200 


2.037(51) 


3.330(225) 


1.770(59) 




0.34510 


0.930(18) 


2.163(75) 


1.764(59) 




34560 


768(121 


1.525(84) 


1.824(51) 




34580 


528 I'M 1 ) 


1.473(102) 


1.902(71) 




34592 


495(141 


1.461(78) 


1.869(51) 




34596 


453(141 


1.722(131) 


1.899(63) 




34600 


o 399(531 


1.542(53) 


1.905(51) 




34604 


423(141 


1.491(35) 


1.851(48) 




0.34608 


0.354(12) 


1.611(35) 


1.809(107) 






0++ 


0++* + gau 


ge d.o.f. 




0.34612 


0.345(6) 


1.710(60) 




0.34616 


0.342(12) 


1.755(87) 




0.34620 


0.318(7) 


1.698(63) 




0.34628 


0.390(6) 


1.656(99) 




0.34640 


0.507(6) 


1.857(108) 




0.34660 


0.567(12) 


2.055(87) 




0.34700 


0.735(9) 


2.181(66) 






0++ 


0++* 


Higgs 


0.37000 


2.723(48) 


5.250(339) 
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Table 14. Masses of 1 and 2 ++ states in units of p| at Mh = 
100 GeV (crossover region) depending on (5h 





Ph 


1~ 


l — » 


-j^ ** 


2++ 


2++* 


symmetric 


0.34200 


3.261(177) 






3.255(84) 






0.34510 


1.374(48) 


2.505(185) 


3.396(252) 


2.121(96) 


2.799(300) 




0.34560 


1.239(30) 


2.304(152) 


3.540(206) 


1.866(63) 


3.057(284) 




0.34580 


1.077(14) 


1.977(102) 


2.919(116) 


1.776(98) 


2.409(201) 




0.34592 


0.948(21) 


1.863(68) 


3.120(162) 


1.670(35) 


2.568(164) 




0.34596 


0.945(21) 


1.854(71) 


3.465(197) 


1.821(53) 


2.733(231) 




0.34600 


0.930(12) 


1.740(135) 


3.258(164) 


1.794(51) 


2.574(147) 




0.34604 


0.816(35) 


1.725(59) 


2.952(137) 


1.692(45) 


2.577(137) 




0.34608 


0.828(26) 


1.581(123) 




1.554(51) 


2.658(192) 




0.34612 


0.738(24) 


1.371(194) 




1.614(29) 


2.592(135) 




0.34616 


0.783(17) 


1.533(132) 




1.566(42) 


2.214(177) 




0.34620 


0.717(32) 


1.833(60) 




1.527(42) 


2.970(210) 




0.34628 


0.684(21) 


2.106(107) 


3.438(275) 


1.608(44) 


2.658(177) 




0.34640 


0.705(26) 


2.301(150) 


3.378(231) 


1.593(56) 


2.433(197) 




0.34660 


0.738(24) 


2.445(63) 


3.402(305) 


1.833(65) 


3.063(300) 




0.34700 


0.867(12) 


2.544(75) 


3.957(357) 


2.013(185) 




Higgs 


0.37000 


2.400(33) 
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